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ABSTRACT 

Recent studies have identified several young stellar objects (YSOs) which exhibit significant mid- 
infrared (mid-IR) variability. A wide range of physical mechanisms may be responsible for these 
variations, including changes in a YSO's accretion rate or in the extinction or emission from the 
inner disk. We have obtained and analyzed multi-epoch near-infrared (NIR) spectra for five actively 
accreting YSOs in the p Oph star-forming region along with contemporaneous mid-IR light curves 
obtained as part of the YSOVAR Spitzer /IRAC survey. Four of the five YSOs exhibit mid-IR light 
curves with modest (~ 0.2-0.4 mag) but statistically significant variations over our 40-day observation 
window. Measuring the strengths of prominent photospheric absorption lines and accretion sensitive 
H I and He I lines in each NIR spectrum, we derive estimates of each YSO's spectral type, effective 
temperature (T e f[), and H band extinction (Ah), and analyze the time evolution of their NIR veiling 
(m and tk) and mass accretion rates (M acc ). Defining a YSO's evolutionary stage such that heavily 
veiled, high accretion rate objects are less evolved than those with lower levels of veiling and ongoing 
accretion, we infer that GY 314 is the most evolved YSO in our sample, with GY 308 and GY 292 at 
progressively earlier evolutionary stages. Leveraging our multi-epoch, multi-wavelength dataset, we 
detect significant variations in mass accretion rates over timescales of days to weeks, but find that 
extinction levels in these YSOs remain relatively constant. We find no correlation between these YSO 
mid-IR light curves and time-resolved veiling or mass accretion rates, such that we are unable to link 
their mid-IR variability with physical processes localized near the inner edge of the circumstellar disk 
or within regions which are directly responsive to mass accretion. We do find, however, that redshiftcd 
He I A10830 emission, where present in our spectra, shows both quantitative and qualitative temporal 
correlations with accretion-sensitive H I emission lines. Blueshifted He I absorption, on the other hand, 
does not demonstrate a similar correlation, although the time-averaged strength of this blueshifted 
absorption is correlated with the time-averaged accretion rate in our sample of YSOs. 
Subject headings: stars: formation - stars: pre-main sequence - stars: variable 



1. INTRODUCTION 

Young stars have long been known to demonstrat e 
significant photometric variability (e.g., Uovl 119451 ). 
Variability at optical and near-infrared (NIR) wave- 
lengths is often attributed to phenomena on or near 
the photosphere, such as hot and cold spots pro- 
duced by accretion and stellar magnetic fields (e.g., 
Herbst et al.lll994tlJohns fc Basri|[l995l: ICarpenter et alj 
20011 : lAlencar et alJ ~ l2001l iBouvier et al.l 120071 ). rT 
cent observations of young stellar objects (YSOs) have 
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also identified si gnificant variabi li ty in the mid-infrared 
(mid-IR; e.g.. iLiu et al.1 119961: IBarsonv et all 120051: 
Morales-Calderon et alj 120091 120101: IMuzer olle et alj 
2009t iFlahertv et al.l 120111 I2012D . where a YSO's spec- 
tral energy distribution (SED) is thought to be domi- 
nated by emission from a circumstellar disk. Dramatic 
changes have been seen over timescales as short as days, 
including v ariations in the shape of the infrared disk 
continuum (lJuhasz et al.1 120071 iMorales-Calderon et al.1 
I2009t IMuzerolle et al.l 120091 : IFlahertv et al.ll2011[ ) and in 
the silicate featur e, which is expected to trace the dust 
in the inner disk (IBarsonv et al.lH997t iSitko et al.l 120081: 
iBarv et al.ll2009t iSkemer et al.ll2010TL 

These recent detections of YSO mid-IR variability have 
been driven by the sensitivity and angu lar resolution 
provi ded by the Spitzer Space Telescope (jWerner et al.l 
12001 . While Spitzer is now in the post-cryogenic phase 
of its mission, the 3.6 and 4.5 ixm channe ls of its InfraRed 
Array Camera (IRAG lFazio et al.ll2004l) continue to pro- 
vide high sensitivity to emission from the inner regions 
of primordial circumstellar disks. The YSOVAR (Young 
Stellar Object VARiability) Spitzer Exploration Science 
Program (PI: Stauffer) is the first large-scale systematic 
survey of YSO variability in the mid-IR. The recently 
published YSOVAR light curves of ~ 2000 objects in 
the Orion Nebula Cluster demonstrate a wide range of 
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morphologies; many YSOs exhibit aperiodic variability 
in the mid-IR, such that evolutionary state appears cor- 
related with both the likelihood and amplitude of vari- 
ability. Class I objects are most likely to be variable, 
and exhibit the largest amplitudes, with Class II YSOs 
(a.k.a. classical T Tauri stars, or CTTSs) and Class III 
YSOs exhibiting progressively fewer , and weaker, varia- 
tions (jMorales-Calderon et"aLll201lD . 

Several scenarios have been suggested to explain 
YSO mid-IR var iability, such as variable mass accretion 
(iVorobvovl 120091 ) . changes in extinction ([Herbst et al.l 
1994) , larg e stru ctural changes in the inner disk 

(|Sitko et al.l 120081 ), changes in the height of the 

disk's inner wall ([Espaillat et al.l 120111 ). stellar oc- 
cultation by density enhancem ents in the inner disk 
( Morales-Calderon et al.1 120111 ). and disk turbulence 
(jTurner et al.1120101 ) . Discriminating between these mod- 
els, however, requires that the mid-IR light curves be 
supplemented with contemporaneous multi-wavelength 
observations. If accretion rate variations drive the ob- 
served mid-IR variability, for example, the strength of 
accretion-sensitive spectral lines should show changes 
commensurate w it h the objects' light curves; indeed, 
iLuhman &: Riekel (|1999l ) (hereafter LR99) observed si- 
multaneous variability in Br7 emission and K band veil- 
ing in two YSOs in p Oph, hinting that accretion activ- 
ity and disk emission may be connected in young stars. 
Extinction and/or emission from localized structures in 
the disk may also contribute to mid-IR variability; simul- 
taneous mid-IR photometry and NIR spectroscopy, such 
as the observations we present in this study, allow these 
scenarios to be tested. 

p Oph is one o f the nearest (d ~ 120 ± 5 pc; 
lLoinard et al.l l2008h star-forming regions, and has been 
extensi vely catalogued in the optical and NIR (e.g. 
LR9 9: iBontemps et all 120011: iWilking et all 120011 : see 
also[Wilking et al. 2008] and references therein). YSOs 
in p Oph are notably variable at NIR wavelengths, 
with many Oph YSOs exhibiting aperiodic variabil- 
ity in the NIR on times c ales of mon t hs to years 
Ce.g. iPlavchan et al.1 l200l iSchobd l20ll iParks et al.1 
lin prepl )": variability at longer wavelengths and on 
shorter timescales, however, is less well documented. In 
this work, we present a coordinated, multi-wavelength 
study of several active YSOs in p Oph to explore physical 
mechanisms that may drive their mid-IR variability. In 
SJ2J we describe our photometric and spectroscopic obser- 
vations and present the resultant light curves and spectra 
for the YSOs in our sample. In $3]we infer spectral types 
and extinction, report time series measurements of each 
YSO's veiling, extinction, and accretion rate, and ana- 
lyze NIR/mid-IR SEDs to study disk emission. Finally, 
we discuss the implications of these results for under- 
standing YSO variability, accretion, and evolution in <21 

2. OBSERVATIONS 
2.1. Mid-IR Photometry 

1 Veiling in the infrared is a measure of excess emission (generally 
interpreted as being due to the presence of a circumstellar disk), 
and is typically defined as r x = -^excess, A /F*, \, where 

^excess, A and 

X are the excess continuum flux (i.e. disk emission) and stellar 
flux at wavelength A, respectively. 



Spitzer /IRAC light curves were obtained for YSOs in 
the p Oph star forming region as part of the YSOVAR 
G06 Exploration Science program. IRAC 3.6 and 4.5 pm 
observations were obtained for three non-contiguous sub- 
fields in the p Oph cluster between 11 April and 17 May 
2010 (MJD 55298-55334); each 5.2' x 5.2' dual-band sub- 
field was centered on an overdensity of previously known 
cluster members. A single observing sequence consisted 
of 8 visits, logarithmically spaced over 3.5 days to sam- 
ple a variety of timescales and to avoid period aliasing 
in characterizing variability; this 3.5 day sequence was 
repeated 12 times to obtain monitoring data spanning 
p Oph's full 40 day Spitzer visibility window. The obser- 
vations were taken in High Dynamic Range mode, with 
a short (0.6 s) and long (10.4 s) exposure taken at each 
position within a five-point gaussian dither pattern. 

We reduced the IRAC images with the custom YSO- 
VAR pipeline, a mo dified version of t h e algo rithm orig- 
inally developed by iGutermuth et al.l (|2009l ) to process 
IRAC data for a large sample of embedded clusters. The 
version of the YSOVAR pipeline used to generate the 
light curves presented here takes as inputs the Basic 
Calibrated Data (BCD) images released by the Spitzer 
Science Center, mosaics the individual images at each 
epoch, and then performs source detection and aper- 
ture photometry on the combined cluster mosaics. The 
sets of long and short exposures at each epoch are sepa- 
rately mosaicked and then combined; however, for bright 
sources with pixels above the nominal saturation limit in 
any long exposure BCD frames, only the short exposure 
mosaics are used for photometry. Sources are then loaded 
into the YSOVAR databas^l, where light curves can be 
constructed for individual objects by merging all detec- 
tions with a 1" search radius, allowing for ^single-pixel 
centroid variations (IRAC pixels are 1.2" in size). The 
resulting IRAC light curves for our spectroscopic targets 
are shown in Figure [1] Typical photometric uncertainties 
for our IRAC measurements are 0.01 to 0.02 magnitudes 
- smaller than the plotted symbols in Figure [TJ To fur- 
ther demonstrate the stability of these observations, we 
include in the figure the light curve of ISY_J162625.24- 
242323.9, a much fainter, non- variable object within the 
p Oph YSOVAR fields. 

2.2. NIR Spectra 

We selected targets for spectroscopic monitoring dur- 
ing the YSOVAR p Oph campaign from the catalog of 
YSOs identified in previous I RAC observations of the re- 
gion (|Gutermuth et al.ll2009D . which identifies 46 YSOs 
in the three YSOVAR/IRAC fields-of-view. Limiting the 
sample to YSOs with J < 11.5 to ensure signal-to-noise 
(S/N) > 50 across the full JHK spectrum in a typical 
5 minute exposure, as required to provide high confi- 
dence detections of relatively weak photospheric absorp- 
tion features and measure veiling from 1-2.5 pm, further 
narrrows the candidates for spectroscopic monitoring to 
five objects. Of these, only the Class II YSOs GY 292 
and GY 314 have shown evidence of active accretion (i.e., 
a positive detection of B r7 or Pa/3 in spect roscopic obser- 
vations by LR99 and/or lNatta et al"1l2006fl : furthermore, 
these two YSOs also demonstrated the potential for ac- 

2 Data for the six objects analyzed in this paper are publicly 
available at http : //ysovar . ipac . caltech . edu/ 
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TABLE 1 

Basic Data for Target YSOs in p Oph 







RA 


Dec 


J 


H 


K S 


Source ID 


other IDs 


(J2000) 


(J2000) 


(2MASS mag.) 


GY 292 


ISO 155 


16:27:33.11 


-24:41:15.3 


11.32 


9.13 


7.81 


GY 308 


IRS 49, ISO 163 


16:27:38.32 


-24:36:58.6 


11.38 


9.43 


8.27 


GY 314 


ISO 166 


16:27:39.43 


-24:39:15.5 


10.75 


9.21 


8.46 


YLW 15a 


IRS 43, ISO 141, GY 265 


16:27:26.94 


-24:40:50.8 


18.53 


13.52 


9.75 


YLW 16b 


GY 274, ISO 145 


16:27:29.41 


-24:39:17.0 


16.72 


12.64 


9.84 



6.6 
6.8 
7.0 
7.2 
7.4 
7.6 

7.8 
8.0 



6.6 
6.8 
7.0 
7.2 
7.4 
7.6 

7.8 
8.0 
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L| C [3.6] 1 
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MJD - 55000 
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to 
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55000 
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310 
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Fig. 1. — IRAC [3.6] (purple circles, left y-axes) and [4.5] (red circles, right y-axes) light curves for the 5 YSOs analyzed here, 
covering the 40-day Spitzer observation window from April 11 to May 17, 2010 (MJD 55298-55334). The left and right y-axes 
are shifted such that the median magnitudes in the two IRAC bands coincide. Typical photometric errors are on the order 
of 0.01 to 0.02 magnitudes, smaller than the plotted symbols in the figure. The three YSOs in the left panel show modest 
(~0.2-0.4 mag), but significant, mid-IR variations, while GY 292 and YLW 15a (right panels) are only marginally variable (see 
also Figure |3j). The light curve of ISY_J162625.24-242323.9, a nonvariable star in the YSOVAR p Oph fields which is notably 
fainter than any of our science targets, is also shown to demonstrate the stability of our calibrated YSOVAR data. 

cretion variability as suggested by two epochs of data 
from LR99. In addition, a third Class II YSO, GY 308, 
also met the S/N cut and showed evidence of both ac- 
cretion activity and variability. While GY 308 did not 
fall directly within one of the dual-band YSOVAR fields, 
it was nonetheless observed (though only in the 3.6 /mi 
channel) due to Spitzer 's simultaneous, offset dual-band 
setup. We thus adopted GY 308 as a third primary tar- 
get, as well as adding two Class I YSOs, YLW 15a and 
YLW 16b, from the dual-band fields as secondary targets 
that were observed when time permitted. Basic proper- 
ties of our five YSOs are presented in Table Q] 



We used the SpeX spectrograph (jRavner et al.l F2003) 
at NASA's 3.0 meter Infrared Telescope Facility (IRTF) 
to observe our targets on 10 nights between 27 March 
and 11 May 2010 (MJD 55283-55328). We present a log 
of our observations in Tabled typical seeing was 
but conditions were rarely photometric, with cloud cover 
ranging from clear to heavily overcast. We conducted our 
observations with SpeX's SXD mode, providing nearly 
contiguous spectral coverage from 0.8 to 2.5 microns. We 
typically observed with a 0.5" slit, providing a nominal 
spectral resolution of R~1200 at the center of each order. 
On two nights of particularly good and bad seeing, we 
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Fig. 2. — J and H band (left) and K band (right) spectra of GY 292, 308, and 314, with an artificially reddened spectral standard shown at the bottom of each panel 
for comparison. Spectra are normalized to unity at 1.5 and 2.1 /im for the J/H and K bands, respectively. Vertical lines indicate the rest wavelengths of key features 
from which we infer spectral types, veiling, and mass accretion rates. Absorption lines appear weaker in the target spectra than in the standards, revealing the presence 
of substantial veiling flux in the target spectra; the difference in slope between the standard and target spectra is primarily due to the fact that the level of veiling is 
wavelength-dependent (see also Figure |4}. Temperature-sensitive photospheric absorption lines show little evidence for variability between epochs; accretion-sensitive 
H I emission lines, by contrast, do appear to change significantly from epoch to epoch. 
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instead utilized the 0.3" and 0.7" slits, providing nominal 
spectral resolutions of R~2000 and ^850, respectively. 

All data were reduced usin g SpeX's dedicated I DL re- 
duction package, SpeXtool (jCushing et al.1 I2004I ) . Tel- 
luric absorption was removed using observations of 
AOV stars at airmasses similar to those of the target 
(A(sec z) < 0.1). Hydrogen features were removed from 
the AOV spectra by dividing by a model of Vega, and 
the target's true spectral slope and absolute flux densi- 
ties were recovered by multiplying the target spectrum 
by a ^10,000 K blackbo dy spectrum scaled to match that 
of the telluric standard (jVacca et al.ll2003[ ). Table EH lists 
the S/N levels reported by SpeXtool for the J, H and K 
bands of each source's reduced SXD spectrum. We de- 
termined relative flux calibration uncertainties in these 
spectra by comparing independent reductions of a single 
raw spectrum calibrated with different telluric (and thus 
flux) standards. When scaled to the continuum, pairs 
of independent reductions agree to within 1-2% in each 
band, and to within ~ 5% across the entire JHK spec- 
trum. 

Our normalized SpeX NIR spectra for GY 292, 308, 
and 314 are shown in Figure [3] along with an artificially 
reddened spectral standard of the appropriate spectral 
type for each YSO (see ^3.2.11 for our spectral type es- 
timates). The difference in continuum slope between 
the reddened standard and object spectra is due to veil- 
ing flux, which increases towards longer wavelengths (see 
Figured]), in the target spectra. None of the YSO targets 
displayed a significant change in spectral slope over the 
course of the monitoring, suggesting that each YSO's 
veiling remained relatively constant over the course of 
this campaign (see ^3.2.2l for a quantitative measurement 
of veiling and its temporal and wavelength dependence) . 

The secondary targets (YLW 15a and 16b) possess suf- 
ficiently high extinctions that obtaining J and H band 
spectra at useful S/N was prohibitively expensive. More- 
over, we were only able to obtain K band spectra for 
these YSOs on five and four occasions, respectively, with 
only one epoch for each secondary target lying within the 
temporal coverage of the IRAC monitoring campaign. As 
a result, these NIR spectra provide no information con- 
cerning the nature of the mid-IR variability exhibited 
by both sources; for completeness we therefore present 
K band spectra for YLW 15a and YLW 16b in Figure [TT1 
of the appendix, and in the spectral archive associated 
with this paper, but we do not analyze these spectra in 
any detail. 

3. ANALYSIS 

3.1. Mid-infrared Variability 

To quantify the amplitude of mid-IR variability in each 
source over our observational baseline, we present in Fig- 
ure [3] cumulative histograms of the IRAC magnitudes in 
each light curve. Stable sources have steep slopes in this 
figure, and variable sources have shallower slopes, as their 
light curves span a larger range of magnitudes. We define 
the variability amplitude indices A[3.6] and A[4.5] as the 
range in magnitudes spanned between the 5% and 95% 
levels in the cumulative histogram for the appropriate 
IRAC band, i.e. the variability amplitude after discard- 
ing the upper and lower 5% of the data. Excluding the 
extrema from this index guards against overestimating 
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Fig. 3. — Amplitude of IRAC variability: cumulative frac- 
tion of epochs as a function of Amag, where Amag at a given 
epoch is defined as the difference between a YSO's magni- 
tude at that epoch and at the epoch of minimum bright- 
ness. YLW 16b exhibits the largest amplitude mid-IR vari- 
ability (,>0.4 mag), while GY 292 is only marginally variable 
(<,(). 1 mag); the other YSOs demonstrate more modest, but 
still significant, mid-IR variations of 0.1-0.3 mag. These am- 
plitudes (which are listed on the figure as A[3.6] and A[4.5]) 
were computed after discarding the upper and lower 5% of 
the data. 

the variability of stable sources due to individual deviant 
data points. We present A[3.6] and A[4.5] for each of our 
five targets in Tabled] 

YLW 16b demonstrates the largest mid-IR variability 
(A[3.6] = 0.44 mag) of the five YSOs in our sample, 
while GY 292 shows the smallest (A[3.6] = 0.07 mag, 
although GY 292 appears to exhibit a short-duration 
flare event near MJD 55325). The mid-IR variabil- 
ity detected from YLW 16b, GY 314, GY 308, and 
YLW 15a represent 5er excursions, with all but YLW 
15a exhibiting A[3.6] > 0.2 mag. GY 292 is not en- 
tirely quiescent, however, as evidenced by comparing its 
light curve with the photometric stability demonstrated 
by the significantly fainter star from the same YSO- 
VAR fields, IS.J162625. 24-242323.9, which has A[3.6] = 
0.05 mag (see Figures [T] and [3]). Although this small 
data set precludes a more quantitative statistical analy- 
sis, t^ejsej^sultsy^ the statistics inferred 
bv lMorales -Caldc ron et al.1 ()2011l ) from the much larger 
Orion YSOVAR dataset. They found that the majority 
of Class I and II YSOs are variable, with a large fraction 
having amplitudes greater than ~ 0.2 mag. In addition, 
the light curves of the most variable objects in our sam- 
ple, GY 314 and YLW 16b, demonstrate strong morpho- 
logical similarities to those of the irregular variables from 
that study. 

Four of our five targets (all but GY 308) have both 
[3.6] and [4.5] photometry; for these YSOs, we examined 
the changes in the mid-IR color [3.6] — [4.5] over time for 
evidence of a color term, which would suggest variability 
in the line-of-sight extinction during our monitoring pro- 
gram. We define S12 as the observed change in [3.6] — [4.5] 
from the median color, then compare the time evolution 
of this quantity with the predictions of two models: (1) 
the case of "gray" variability, i.e. 812 = 0, and (2) the ex- 
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TABLE 2 

IRTF/SpeX Observing Log and spectrally derived colors 





1 Observation 














Date 


exp. 


SXD 


J 


H 


K 


Star 


(UT, 2010) 


time (s) 


slit (") 


S/N a 


S/N a 


S/N a J - H b H - K s h 



GY 292 


March 27 


360 


0.5 


46 


168 


295 


2.23 


1.26 




April 01 


360 


0.5 


111 


312 


384 


2.19 


1.20 




April 02 


720 


0.5 


162 


452 


603 


2.20 


1.20 




April 09 


720 


0.5 


115 


344 


429 


2.23 


1.26 




April lb 


840 


0.5 


60 


172 


209 


2.16 


1.24 




April 19 


1440 


0.5 




83 


155 


1.93 


1.26 




May 02 


480 


0.5 


75 


261 


382 


2.16 


1.24 




May 03 


480 


0.7 


132 


375 


512 


2.17 


1.28 




May 07 


360 


0.5 


133 


356 


449 


2.18 


1.27 




May 11 


720 


0.3 


132 


334 


391 


2.18 


1.27 


GY 308 


March 27 


360 


0.5 


52 


165 


233 


2.04 


1.07 




April 01 


360 


0.5 


83 


219 


277 


1.99 


1.03 




April 02 


720 


0.5 


148 


366 


446 


1.98 


1.02 




April 09 


720 


0.5 


86 


249 


314 


1.97 


1.00 




April 16 


1080 


0.5 


70 


183 


216 


1.98 


1.05 




April 19 


1080 


0.5 


40 


127 


169 


1.97 


1.01 




May 02 


840 


0.5 


49 


154 


200 


2.00 


0.91 




May 03 


480 


0.7 


123 


312 


413 


1.91 


1.01 




May 07 


360 


0.5 


116 


286 


355 


1.91 


1.00 




May 11 


360 


0.3 


83 


215 


267 


1.91 


1.00 


GY 314 


March 27 


360 


0.5 


82 


184 


224 


1.63 


0.78 




April 01 


720 


0.5 


166 


310 


324 


1.64 


0.79 




April 02 


720 


0.5 


242 


456 


522 


1.60 


0.76 




April 09 


720 


0.5 


170 


333 


375 


1.62 


0.77 




April 16 


1260 


0.5 


116 


280 


367 


1.58 


0.80 




April 19 


1320 


0.5 


124 


227 


227 


1.66 


0.89 




May 02 


960 


0.5 


29 


90 


130 


1.60 


0.78 




May 03 


480 


0.7 


177 


345 


411 


1.57 


0.77 




May 07 


360 


0.5 


173 


323 


365 


1.56 


0.79 




May 11 


720 


0.3 


179 


314 


325 


1.59 


0.81 


YLW 15a 


March 27 


960 


0.5 




5 


158 


3.55 


3.72 




April 01 


1440 


0.5 




5 


165 


5.77 


3.74 




April 02 


1320 


0.5 




11 


276 


5.58 


3.81 




April 09 


1440 


0.5 




9 


240 


4.74 


3.86 




May 07 


600 


0.5 




9 


214 


2.95 


3.75 


YLW 16b 


March 27 


960 


0.5 




87 


311 


2.45 


2.44 




April 01 


1440 


0.5 




71 


361 


2.54 


2.59 




April 02 


960 


0.5 


7 


108 


424 


3.74 


2.57 




May 07 


720 


0.5 


13 


131 


362 


3.49 


2.30 



a Signal- to- noise ratios (S/N) as reported by 
b NIR colors were computed by convolving our 



tinction law of llndebetouw et all (|2005l ). which predicts 
5 12 = {1-A 2 /A 1 )S 1 , where ^ 2 Mi = 0.77±0.16 is the dif- 
ferential extinction in the two IRAC bands derived from 
Table 1 of that work, and <5i is the change in [3.6] from 
the median. Computing the reduced \ 2 goodness-of-fit 
statistic for the data as compared with each model, we 
find the following: for the two YSOs with more modest 
mid-IR variability (GY 292 and YLW 15a), both mod- 
els are compatible with the data (\ 2 < 2); this simply 
reflects the fact that the uncertainties in S12 are similar 
in scale to the level of variability. Conversely, neither 
model is particularly consistent with the observed color 
changes (x 2 > 6) for YLW 16b, which shows stochastic, 
measurable changes in 612 and no evidence for a color 
term. Finally, for GY 314, the extinction model is more 
consistent with the data (Xcxt = 1-4) than the gray model 
(Xgray = 5-0), but the trend appears to be leveraged by a 
small subset of the data, with the vast majority actually 
clustered near (<5i,<5i 2 ) = (0,0). Given that xlxt > 1; 
the relatively high uncertainty on the extinction model 
slope, and the fact that the entire range in 612 is less than 
0.2 mag for GY 314, we are unable to confidently claim 



SpeXtool. 

SpeX spectra with standard 2MASS passbands. 

that this extinction law accurately reproduces the data 
for this YSO. We thus see no evidence for (or against) 
extinction-induced mid-IR variability in any of our tar- 
gets from their light curves alone; in Ej4.2.2[ we utilize 
our spectral data to investigate this possibility in more 
detail. 



3.2. Spectral Measurements 

A few key parameters dominate the morphology of a 
YSO's NIR spectrum: (1) the effective temperature of 
the stellar photosphere (T c ff), primarily diagnosed by the 
presence and relative strength of prominent (Na, CO, Ca, 
H 2 0, etc.) absorption features; (2) the amount of veil- 
ing flux, diagnosed by the absolute depth of photospheric 
features relative to a template of similar spectral type; 
(3) extinction, diagnosed by the broadband slope of the 
spectral continuum; and, (4) the mass accretion rate, di- 
agnosed by the strength of prominent emission lines such 
as Paschen /3 (Pa/J, H I A12818) and Brackett 7 (Br7, 
H I A21661). In this section, we follow previously em- 
ployed techniques to infer these physical parameters for 
our YSO targets from moderate resolution NIR spectra 
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TABLE 3 

MlD-IR AND NIR VARIABILITY INDICES 



Object 


A [3.6] 
(mag) 


A[4.5] 
(mag) 


(EW Bn ) a 
(-A) 


min. EWfl r7 a 
(-A) 


max. EWs,- 7 a 
(-A) 


(-A) 


{ t h) 


(ric) 


GY 292 


0.07 


0.07 


4.9 ± 1.0 


1.7 ±0.6 


11.6 ±5.1 


5.5 ± 1.1 


1.4 


2.5 


GY 308 


0.18 




3.7 ±0.6 


1.4 ±0.5 


6.4 ± 1.6 


4.4 ±0.7 


0.9 


1.2 


GY 314 


0.21 


0.17 


1.6 ± 0.2 


0.8 ±0.8 


2.5 ±0.8 


0.7 ±0.3 


0.5 


0.9 


YLW 15a 


0.12 


0.09 














YLW 16b 


0.44 


0.42 















a Veiling-corrected EWs 

(e.g., LR99: IWirrston et al.|[2Q0l iCovev etall [20Toh to 

assess which, if any, physical properties correlate strongly 
with mid-IR variability. 

3.2.1. Spectral Types 

Visual inspection of strong, temperature sensitive ab- 
sorption lines in Figure [2] does not reveal any evidence 
for significant changes in the photospheric temperature 
of these YSOs. Specifically, the relative strengths of 
nearby temperature sensitive lines (e.g., the 2.21 /im 
Na I & 2.26 fim Ca I doublets) appear consistent from 
epoch to epoch, modulo small changes in overall depth 
due to modest changes in the star's continuum veiling. 
Given the evident stability in each YSO's photospheric 
temperature, we have estimated each star's photospheric 
spectral type from the spectrum obtained on April 1st; 
we consider potential veiling and extinction variations 
in more detail in §§ 13.2.21 and 14.2.21 respectively. We 
estimated each star's spectral type using the algorithm 
developed and presented in full by ICovev et al.l (|2010f ): 
we refer the interested reader to that work for a com- 
plete description of the spectral analysis procedure, and 
simply summarize here the most pertinent details. 

Spectral type and veiling estimates are derived by 
comparing the absolute and relative strengths of tem- 
perature sensitive spectral features (e.g., Na, Ca, CO, 
Mg, H 2 Q; for a full list see Tables 3 and 4 of 
ICovev et al.ll2010|) measured from spectral standa r ds ob - 
serv ed bvlLancon fc Woodl (|2000f ). iCushing et al.l (|2005l ) 
and ICovev et al.l (|2010l )~ These standards span a wide 
range of effective temperatures, and include giants and 
dwarfs, as well as pre-main sequence stars in the TW 
Hya association, enabling an estimation of the target's 
surface gravity as well as its T c g. Estimates of the ex- 
tinction toward each YSO photosphere are determined 
by comparing its NIR spectral slope, as quantified by its 
spectrally integrated J — H and H — K s colors, to the 
colors expected for an artificially reddened and veiled 
spectral standard. We computed these spectrally inte- 
grated colors by convolving each SpeX spectrum with 
standard 2MASS filter curves and integrating the total 
transmitted flux, u sing routines originally developed by 
ICovev et al.l (|2007[ ) . The spectrally integrated J — H and 
H — K s colors derived from each target spectrum are re- 
ported in Table El 

We determined the spectral type, veiling and extinc- 
tion that provide the best overall fit to each YSO's NIR 
spectrum using an iterative technique to refine the ini- 
tial values we assume for each parameter. We began by 
adopting a fiducial K7 spectral type for each YSO, a com- 
mon spectral type for typical T Tauri stars. Comparing 



the overall spectral slope (as parameterized by the spec- 
trally integrated J — K s color) and strengths of promi- 
nent absorption features in the target spectrum to those 
of a pre-main sequence K7 spectral standard, we derived 
initial estimates of the YSO's extinction and veiling, re- 
spectively. We then synthetically reddened and veiled 
the grid of spectral standards to match these initial es- 
timates, and compared the strengths of the temperature 
and gravity-sensitive absorption features in the YSO's 
spectrum to those measured from the synthetically red- 
dened and veiled standards. The un-reddened, un-veiled 
spectrum of the new T e g standard could then be used 
to produce updated estimates of the source's veiling and 
extinction, which in turn could be re-applied to the grid 
of standard spectra. Repeating this process several times 
was sufficient to identify a self-consistent set of extinction 
(which we report in the H band as Ah), veiling, and spec- 
tral type estimates that reproduced each targets spectral 
morphology on both large and small spectral scales. We 
compute T e g from spectral type using the scale derived 
by LR99, and report the spectral type (SpT), T e g, and 
Ah for each of our targets in Table |6] alongside our de- 
rived mass accretion rates (see £13.2.31) . 

Previous tests of this algorithm have identified a char- 
acteristic uncertainty in spectral type estimates of <r=±l 
subclass; to incorporate the possibility of veiling and/or 
extinction variability, we conservatively adopt an uncer- 
tainty of (T=±2 subclasses in our spectral type deter- 
minations. Estimating spectral types for YSOs, partic- 
ularly at the youngest evolutionary stages, is challeng- 
ing, and there are discrepancies between the previously 
published results for these YSOs. LR99 used moderate- 
resolution K band spectra to derive spectral types of K8, 
K8, and MO for GY 292, 308, and 314, respectively, while 
iGatti et all (|2006[ ) determined somewhat earlier values 
of K3, K5, and K5 from low-resolution J band spectra; 
both authors report uncertainties of ^3 subclasses. Our 
results of K5, M0, and M2, respectively, agree with the 
types LR99 derived for all three YSOs to within our ±2 
subclass uncertainty, but a re in somewhat wo rse agree- 
ment with those derived by IGatti et al.l (|2006l) . 

We suspect that the agreement between the types de- 
rived here and by LR99, and the s ystematic off s et to - 
wards earlier types in the analysis of IGatti et al.1 (2006) 
is primarily due to a key difference in the treatment of 
veiling in these works. The types inferred here and by 
LR99 explic i tly all ow for line dilution by veiling, whereas 
IGatti et all {2006) assume that the J band is veiling- 
free, and derive types based on absolute line strengths. 
If these YSOs possess non-ne gligbl e J band veiling, as 
many CTTSs do (jCieza et al.l l2005D . their J band pho- 
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tospheric lines will appear systematically weaker than 
they would in the absence of veiling. Many atomic fea- 
tu res in NIR s tellar spectra, including those analyzed 
bv lGatti et al.1 (|2006| ). weaken towards hotter T e gs, such 
that the dilution of a YSO's line strengths by veiling 
flux can be misdiagnosed as indicative of a warmer pho- 
tosphere, and thus an earlier spectral type. The analy- 
sis employed here, as well as by LR99, assigns spectral 
types based on ratios of lines with somewhat divergent 
temperature dependences, and thus minimizes (but does 
not completely remove) this degeneracy between T e g and 
veiling flux; once the YSO's spectral type has been estab- 
lished by these line ratios, the absolute line depths can 
be used to infer the presence of additional veiling flux. 
Furthermore, the trend between absolute line strength 
an d T s ff is rather sha llow for the J band features used 
bv iGatti et al.l (|2006l ) to infer spectral types, such that 
the intrinsic scatter between line strengths at a particular 
T e ff is comparable to the m agnitude of the slow trend over 
a full spectral class. If the IGatti et al.l (|2006f ) uncertain- 
ties were increased to reflect this additional degeneracy, 
it would better reconcile their results with ours. 

The degeneracy between T e g and veiling discussed 
above is sufficiently challenging that veiling remains the 
most poorly constrained parameter by our iterative spec- 
tral fitting technique. Reasonable by-eye fits to the ob- 
ject spectra can be obtained with veilings that differ by 
as much as 50%, after allowing for a slight modification 
in T e ff (which we incorporate into our spectral type un- 
certainty) to offset the change in absolute line depth. We 
therefore describe in the next section a separate, objec- 
tive measurement of potential veiling variations in each 
YSO's spectrum once an appropriate spectral template 
has been identified. 

We estimate the uncertainty in the extinction An 
derived from our spectral analysis by deriving alter- 
native An estimates from dereddening the published 
2M ASS magnitudes f or our YSOs to the CTTS lo- 
cus (jMever et al.lll997ft . Comparing these 2MASS- based 
An values to those obtained using our spectrally inte- 
grated JHK S magnitudes reveals that the two estimates 
are consistent to ^0.3 magnitudes in An, and the dif- 
ferences may be partially due to long-term photomet- 
ric variability. However, as further discussed in i)3.31 
the NIR magnitudes of these YSOs are r elatively stable 
betwe en the measurements obtained by iBarsony et al.l 
(|1997t ) and the 2MASS survey, so we conservatively adopt 
the full 0.3 mag as the uncertainty in An ■ For consistency 
with the line luminosity calculations in section ["3.2.3[ the 
Ah values listed in Table |5] are those determined from 
their published 2MASS magnitudes. 

3.2.2. Veiling Measurements 

As mentioned in the previous sections, young stars of- 
ten exhibit 'veiling', an excess of near- infrared emission 
which presumably arises from the inner edge of a cool 
circumstellar disk. As a result, photospheric absorption 
lines in YSO spectra can be significantly weaker than in 
field stars of similar T e g when the spectra are normal- 
ized to the continuum. In this section, we analyze our 
multi-epoch NIR spectra to directly measure the time 
evolution of NIR veiling in our spectroscopic targets. 

To begin, we parameterize veiling in terms of an ab- 
sorption feature's equivalent width (EW), a continuum- 



independent observable that is robust against absolute 
flux calibration uncertainties. If we define In nc to be the 
total absorbed flux in a photospheric line, the EW of a 
spectral feature in an unveiled star's spectrum is given by 
EW ull vciicd = Iiinc/(F*) , where (F*) is the continuum flux 
at the wavelength of the feature. In contrast, the contin- 
uum flux in a star with the same T e ff but nonzero excess 
emission from the disk is F cont — F cxccss + i 7 "*, where 
^excess is the disk emission at that wavelength. Thus the 
EW of the exact same spectral line in the veiled star will 
instead be EW vci i od = hne/ (F*+F excess ) . A target star's 
veiling r\ can therefore be calculated from the equivalent 
widths of an absorption line in the spectrum of the target 
and in an unveiled spectral standard of the same spectral 
type. Using the above relations with the fact that veiling 
is defined as r\ = F excess /F* yields the relation: 



To estimate each target's veiling at all epochs and a range 
of wavelengths, we measured the EWs of several promi- 
nent H and K band photospheric absorption features 
in each spectrum, selecting consistent line centers and 
widths for each feature so as to avoid contamination from 
neighboring features while still capturing the full line pro- 
file at each epoch. We also chose the sizes and widths 
of continuum regions on either side of each line to avoid 
nearby lines. We present the complete list of features 
measured in each target (including the emission features 
discussed in ^3.2.31) in Table [4] along with our chosen line 
and continuum parameters. We also measured the EWs 
of each absorption feature in standard stars of the ap- 
propriate spectral type, then used Equation [T] to infer 
veiling for each of our targets at each epoch. Table [3] 
presents the minimum, maximum and mean veilings, as 
well as the dispersion in the veiling time series, for each 
measured feature^. 

Emission from a circumstellar disk typically increases 
or remain s relatively constant th rough the NIR and mid- 
IR0 (e.g., iD'Alessio etaLI I2006D . in contrast to photo- 
spheric emission which is decreasing along the power- 
law Rayleigh- Jeans blackbody tail at these wavelengths. 
Thus one would expect veiling in YSOs to increase across 

3 Our measured K band veilings are systematically higher 
than those of LR99 by a factor of ~ 2. To investigate this dis- 
crepancy, we compared veiling measurements obtained using dif- 
ferent combinations of target spectra and measurement methodol- 
ogy. Specifically, we compared veiling measurements obtained by: 
(1) applying the methodology described above to the LR99 spec- 
tra, and (2) applying the LR99 methodology to our SpeX spec- 
tra (K. Luhman, private communication). The veilings we derive 
from these tests are consistent with those in Table [5] to well within 
our uncertainties. This suggests the most likely cause for the dis- 
agreement between our measurements and those of LR99 lies in 
the choice of spectral standards used to define an unveiled photo- 
sphere: the veilings listed in Table [5] and via these tests were all 
derived using SpeX spectra of dwarf stars as proxies for unveiled 
photospheres, whereas LR99 measured veiling with respect to spec- 
tral standards generated by averaging spectral of dwarf and giant 
stars. 

4 Disk SEDs are in truth exceptionally complicated in the in- 
frared; the exact details depend on a host of factors including 
dust grain size distribution, grain growth properties, disk geom- 
etry, viewing angle, stellar irradiation, and interstellar UV field. 
Our goal is only to empirically measure the wavelength dependence 
of veiling in our targets without attempting to model the disk in 
any great detail. 



Oph SpeX Variability 



9 



TABLE 4 
Spectral Line Parameters 



Line 


Line Center 


Line Width 


Cont. l a 


Cont. l a 


Cont. 2 a 


Cont. 2 a 


Name 


(unci) 


(am) 


Center (/mi) 


Width (um) 


Center (|im) 


Width (/im) 


Mg I 


1.5037 b 


0.005 


1.4982 


0.0028 


1.5091 


0.0028 


Al I 


1.6740 


0.006 


1.6696 


0.0020 


1.6781 


0.0020 


Mg I 


1.7110 


0.003 


1.7044 


0.0032 


1.7177 


0.0032 


Ca I 


1.9820 c 


0.013 


1.9681 


0.0027 


1.9972 


0.0027 


Na I 


2.2071 b 


0.006 


2.1991 


0.0064 


2.2152 


0.0064 


Ca I 


2.2640 b 


0.008 


2.2565 


0.0043 


2.2715 


0.0043 


Mg I 


2.2812 


0.003 


2.2755 


0.0027 


2.2868 


0.0027 


He Wuc d 


1.08220 


0.0014 


1.0803 


0.0013 


1.0859 


0.0013 


He re d e 


1.08335 


0.0009 


1.0800 


0.0013 


1.0859 


0.0013 


Pa7 


1.0940 


0.0025 


1.0908 


0.0019 


1.0972 


0.0019 


Pa/3 


1.2817 


0.0025 


1.2782 


0.0023 


1.2853 


0.0023 


Br7 


2.1655 


0.006 


2.1586 


0.0038 


2.1725 


0.0038 



1 Continuum regions were chosen so as to include the entirety of the spectral line at all epochs for all 
YSOs while minimizing contamination from neighboring features. 

2 Doublet; the reported line center is the mean of the two individual line centers 

3 Triplet; the reported line center is the mean of the three individual line centers! 

4 The portion of the He line just blueward of the star's rest frame velocity (see i|4.3H 

5 The portion of the He line just redward of the star's rest frame velocity (sec i| 1.311 



TABLE 5 

Veiling measured from photospheric absorption lines 









GY 292 






GY 308 






GY 314 




line 


\(fim) 


^*A,min 


f'A,max 




(ay) 


' A, min 


"A, max 


(rx) 


(Or) 


^*A,min 


^ A,max 




(oy) 


Mgl 


1.5037 


0.58 


1.00 


0.81 


0.15 


0.26 


0.54 


0.36 


0.09 


0.00 


0.16 


0.02 


0.08 


Al I 


1.6700 


0.96 


1.84 


1.46 


0.32 


0.40 


1.49 


0.95 


0.35 


0.42 


1.10 


0.68 


0.26 


Mgl 


1.7110 


0.87 


2.39 


1.40 


0.43 


0.66 


1.11 


0.86 


0.18 


0.33 


0.73 


0.52 


0.14 


Ca I 


1.9820 


1.13 


3.23 


1.99 


0.71 


0.78 


1.58 


1.05 


0.23 


0.56 


1.40 


0.97 


0.31 


Na I 


2.2071 


1.64 


3.44 


2.49 


0.68 


0.80 


1.64 


1.22 


0.27 


0.69 


1.36 


1.01 


0.24 


Ca I 


2.2640 


1.58 


3.81 


2.54 


0.80 


0.80 


1.50 


1.14 


0.24 


0.59 


1.18 


0.82 


0.22 


Mgl 


2.2820 


0.90 


3.80 


1.89 


0.95 


0.56 


1.53 


0.91 


0.32 


0.27 


0.62 


0.41 


0.11 



the NIR. To confirm this scenario, we computed time- 
averaged line-by-line veiling measurements for each YSO, 
which are presented as a function of wavelength in Fig- 
ure IH solid error bars show average measurement errors, 
while dotted error bars represent the formal standard 
deviation of the veiling time series. Veiling does indeed 
increases with wavelength across the J UK bands in all 
three YSOs, although there is some evidence for a de- 
crease in slope around 2.2 /im. We note, however, that 
the Mg I A22814 feature was particularly weak in our 
targets, leading to high formal uncertainties on veiling 
measurements at 2.28 /im, and thus we do not accord 
any physical significance to the apparent systematic de- 
crease in veiling at the upper end of the K band. 

To characterize the typical veiling levels in the H 
and K bands and their potential time variability for 
comparison with our IRAC light curves, we select the 
strongest features near the center of each band (and 
therefore best separated from telluric H2O features), and 
define the veiling inferred from these lines as rjj and 
tk- In particular, we define m — r(Mg I A17110) and 
rjf = [r(Na I A22071) + r(Ca I A22640)]/2. We present 
time-averaged veilings r# and tk in Tabled and visually 
show the time evolution of tk alongside the other NIR 
and mid-IR indices in Figure[7] Veiling in all three YSOs 
appears to vary only modestly over our observation win- 
dow (note that variations in tk have been suppressed by 
a factor of two in the figure for readability). 



We estimated uncertainties in EW measurements 
through a Monte Carlo approach devised for SpeX NIR 
spectra (see lMuirhead et aLll2011l for a complete descrip- 
tion). In summary, artificial noise was added to each 
YSO spectrum, with the noise drawn from a Gaussian 
probability distribution function (PDF) having a width 
equal to the flux uncertainty estimated for each pixel by 
the SpeXtool pipeline. The EW of each spectral feature 
was then measured from this artificially noised spectrum. 
This process was independently repeated 1000 times for 
each EW measurement; the error in a given EW mea- 
surement is then taken to be the standard deviation of 
all the simulated EW measurements for that feature. To 
independently estimate the effect of a potentially mis- 
placed continuum, we separately computed EWs of sev- 
eral emission and absorption features in our fiducial April 
1 spectrum while artificially varying the continuum level 
by ±2% (the approximate range over which the by-eye 
placement appeared reasonable). Over this range of con- 
tinuum placement levels, EWs varied by an average of 
- 0.10 A, so we added this value in quadrature with 
the Monte Carlo-determined error to determine the total 
EW uncertainty, which is then formally propagated into 
our veiling and mass accretion rate errors. 

We note that in all three YSOs, the photospheric fea- 
tures seen in each spectrum appear to be systematically 
weaker in the last six epochs as compared with the first 
four, leading to relatively higher veiling measurements 
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Fig. 4. — Veiling as a function of wavelength, as inferred from 
individual spectral features in GY 292, 308, and 314. The 
solid error bars show the average measurement uncertainties, 
while the dotted error bars represent the la variability in veil- 
ing across all epochs. Veiling increases as a function of wave- 
length across the J and H bands in each YSO's spectrum, 
then flattens across the K band. The ordering of objects by 
veiling strength is identical across the NIR, with GY 292 and 
314 exhibiting the largest and smallest values, respectively. 
The Mg I A22814 feature is very weak in our targets, and 
thus the uncertainties in veiling inferred from this line are 
large. 

at the end of the monitoring program. This could po- 
tentially be due to the fact that different standard stars 
were used for telluric corrections in the first four epochs 
as compared to the last six. To explore this possibility, we 
consider the relative strengths of absorption features in 
target spectra reduced with both telluric standards, with 
all spectra obtained at comparable airmasses on April 16 
(MJD 2455303). Specifically, we measure the EW of the 
full suite of absorption features listed in Table 0] in each 
of two alternate reductions of both GY 292 and GY 314, 
one reduced with HD 145188 (the standard used for the 
majority of the reductions for the last six epochs), and 
the other reduced with HD 145127 (one of the standards 
used for the first four epochs). We find that EWs mea- 
sured from spectra reduced with HD 145188 are system- 
atically lower (in all lines and both YSOs) by ^0.1-0.5 A 
(l-2tr) than those reduced with HD 145127. As Equa- 
tion[T]shows, a 0.1-0.2 A change in the EW of a feature in 
the standard star's spectrum (EW unV eiied) leads to a 15- 
30% difference in measured veiling in these YSOs. There 
are no obviously visible photospheric absorption features 
in the standard star spectra, but we are unable to rule 
them out below the ~0.5 A level due to the large num- 
ber of nearby telluric features. There is also a substantive 
difference in the characteristic airmasses of the observa- 
tions obtained during the first four epochs (secz = 1.4 
to 1.5) compared to the last six epochs (secz = 1.6 to 
2.2), so imperfect subtraction of atmospheric lines at 
high airmasses could have also contributed to this ef- 
fect. Nonetheless, the IRAC observation window only 
overlaps the last six spectroscopic epochs, and thus our 
primary conclusions (which are based on relative, direct 
temporal comparison of mid-IR light curves and veiling) 



are entirely unaffected by this potential systematic effect. 
In addition, hydrogen features (e.g. Br7 and Pa/3) are 
well-subtracted in the telluric calibration, so EW mea- 
surements of these lines are insensitive to the choice of 
standard star. 

3.2.3. Accretion & Outflow Line Diagnostics 

Muz erolle et al.l ()1998D demonstrated that the lumi- 
nosities of NIR hydrogen emission lines (specifically, Br7 
and Pa/3) emitted by young, low-mass (0.2-0.8 M & ) stars 
correlate well with their mass accretion rates as derived 
from ultraviolet (UV) observations. These authors also 
show a slightly weaker correlation between the equivalent 
widths (EW) of these lines and YSO mass accretion rates; 
the advantage of using EWs for our analysis is that they 
reflect the intrinsic strength of the line relative to the 
local continuum, and thus avoid errors associated with a 
spectrum's absolute flux calibration. As our spectra are 
subject to substantial uncertainties in their absolute flux 
calibration, Br7 and Pa/3 EWs provide the most reliable 
means of investigating changes in our YSO targets' ac- 
cretion rates over time; due to the superior S/N achieved 
in our K band spectra, we selected Br7 in particular for 
detailed time series analysis of the mass accretion rate. 

Our measurements of EW(Br7) at each epoch are 
presented visually alongside our other NIR and mid-IR 
quantities in Figure [7] variability in this feature for all 
three YSOs is clearly evident (note that variations have 
been suppressed by a factor of two in the figure for read- 
ability) . Combining our measurements to study the char- 
acteristic accretion activity of our YSO targets over the 
entire 40-day observation window with improved statis- 
tical power, we also computed time-averaged equivalent 
widths (EW(Br7)) for each YSO, which we present in Ta- 
ble |31 Both individual and time-averaged measurements 
have been veiling-corrected to remove the disk contribu- 
tion to the continuum that would otherwise bias EWs 
downward in more heavily veiled spectra. Note that the 
uncertainty reported on (EW(Bi - 7)) is the formal uncer- 
tainty of the mean, defined as er m = a/*/N, where a is 
the standard deviation of the data and N is the number 
of measurements. 

To place each YSO's characteristic mass accretion 
rate on an absolute physical scale, we converted time- 
averaged Br7 EWs into line luminosities by multiplying 
by the continuum flux level inferred from the 2MASS 
K s m agnitude, after corr ecting for distance (d~ 120 ± 
5 pc; lLoinard et al.ll2008l ) and extinc tion (using our de - 
rived Ah and the extinction law of iFitzpatrickl [1999) . 
The total accretion luminosity L acc of each object can 
then be derived from the empirical correlations be- 
tween L acc and Brr line luminosity first presented by 
iMuzerolle et al.l (|1998| ) and later refined bv iNatta et al.l 
( 20061 see their equations 3 & 4). The mass accretion 
rate is then given by M acc — L acc R*/(GM r ). For this 
calculation, we estimated each YSO's radius by interpo- 
lating the T e g estimates inferred from our spectral typ- 
ing procedure onto t he pre-main s e quenc e evolutionary 
tracks computed by iBaraffe et al.l (|1998t ). after adopt- 
ing a characteristic age of 3.5 Myr (the m edian age of 
p Op h YSOs inferred from this model grid; iCovev et al.l 
120 lOf); for each YSO , we adopted the mass reported 
bv INatta et "all (|200ll . GY 292 and GY 308 also ex- 
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hibit strong Pa/3 emission lines (Figure [5]) , and we follow 
the above procedure using time-averaged Pa/3 EW and 
2MASS J magnitude to derive M aC c from the Pa/3 line 
in these two objects as well. We summarize the H I line 
properties and derived accretion rates in Table |5] the 
implied time-averaged mass accretion rates range from 
about 10" 8 - 5 M Q /yr for GY 292, to 1CT 9 - 2 M©/yr for 
GY 314. The combination of uncertainties in Ah (~ 
0.3 mag), distance (~ 4%), and the EW measurement 
(~10-20%) lead to a formal uncertainty of ~ 30% in 
Br7 and Pa/3 line luminosity. Errors in accretion lumi- 
nosity are thus dominated by the scatte r in th e L acc — line 
luminosity relation, which iNatta et all ([2004) place as a 
factor of 2 to 3. For M acc , uncertainties in mass (~ 60%) 
and radius (~ 40%) also contribute, and we conserva- 
tively estimate a total error of a factor of 3 (equivalent 
to ~ 0.5 dex) in M acc . Despite the relatively large uncer- 
tainties, we see a measurable difference in time-averaged 
mass accretion rages amongst our three YSO targets; we 
will discuss these results in the context of YSO evolution 
in j gZD 

For the two YSOs with detected Pa/3 and Pa7 emis- 
sion (GY 292 and GY 308), we can compute Br7/Pa/3 
and Pa7/Pa/3 line ratios to look for consistency with 
Case B recombination theory, which assumes that the 
emitting gas is optically thick to Lyman series photons 
and optically thin to a ll other H I transitions (e.g., 
iStorev fc Hummer! 119951) . Using the extinctions (Ah) 
derived in %T2~Tl we dereddened each YSO' s 2MASS J 
and K band magnitudes using an R v = 3.1 lFitzpatrickl 
(1999) extinction law. Converting these dereddened 
2MASS magnitudes into estimates of each YSO's de- 
reddened continuum flux densities, we converted our Br7, 
Pa/3, and Pa7 EWs into line fluxes. In Figure [5] we plot 
line fluxes F^ T1 and -Fp a7 as functions of -Fp a( 3 for both 
GY 292 and GY 308. Computing the weighted aver- 
age ratios from the data, we find Br7/Pa/3 line ratios of 
0.30 ± 0.07 and 0.29 ± 0.10, and Pa7/Pa/3 line ratios of 
0.61±0.14 and 0.69±0.35 for GY 292 and GY 308, respec- 
tively. Line ratios are consistent between the two YSOs, 
and there is no evidence for significant variability in ei- 
ther line ratio to within our measurement uncertainties. 
Note that the la errors listed above represent the stan- 
dard deviations of the line ratio time series, and include 
neither measurement uncertainties nor the additional er- 
ror introduced by assuming that the 2MASS magnitudes 
represent accurate continuum fluxes. The uncertainties 
in the calculated Pa7/Pa/3 line ratios are likely to be 
even higher due to the fact that we use a single value 
for the J band continuum (the 2MASS measurement) 
to compute both line fluxes despite their relatively large 
wavelength separation. Our derived Br7/Pa/3 results are 
consistent with the line ratio of ~ 0.25 measu red in a 
sample of YSOs in Taurus bv lBarv et al.l (|2008| ). as well 
as with the predictions of Case B recombination for a 
wide range of realistic electro n temperatures and densi- 
ties (|Storev k, Hummerll 19951 ). However, the above com- 
plications, combined with the fact that our S/N is inade- 
quate to reliably measure additional higher-rt transitions 
in either the Paschen or Brackett series, render us unable 
to directly constrain the electron temperature or density 
in the emitting regions of these YSOs using our observa- 
tions. 



The He I A10830 line has also proven to be an infor- 
mative pro be of outflow and accretion activity in y oung 
stars (e.g. JEdwards et al.ll2006t iFischer et a!T[2008l) . To 
ensure the reliability of the kinematic content of the rich 
He I line profile, we have compensated for any potential 
remaining wavelength calibration errors in our spectra 
by adjusting the Pa/3 line center to the Pa/3 rest wave- 
length at all epochs using the following procedure. First, 
we measured the center of the Pa/3 feature in the spec- 
tra of GY 292 and GY 308 (we excluded GY 314 from 
this procedure because it did not exhibit Pa/3 emission at 
all epochs). We then averaged the offsets measured for 
GY 292 and GY 308 at each epoch, and shifted all spec- 
tra by that amount to center Pa/3 at its rest wavelength. 
The typical shift applied was about 2 A, or roughly half 
a pixel, comparable to the uncertainty in the measured 
Pa/3 line centers. The average difference between the off- 
sets derived for GY 292 and GY 308 was 1.7 A, with no 
systematic behavior in the size or direction of observed 
shifts. Pa/3 is separated from the He A10830 line by a few 
hundred pixels in our spectra, so we also measured the 
offset of the Pa7 feature from its rest wavelength at each 
epoch in order to test the appropriateness of this wave- 
length recalibration (Pa7 is much closer to the He I line, 
but is considerably weaker than Pa/3, providing far less 
accurate wavelength centroids). The difference between 
the Pa/3 and Pa7 offsets was on average about 1.4 A - 
less than the uncertainty on the Pa/3 centroids - demon- 
strating that these wavelength corrections are robust at 
the sub-pixel scale in the wavelength regime of He I. 

Having verified the fidelity of the wavelength solution 
near the He I A10830 line, we characterized the struc- 
ture of the He I line profile in each of our target spectra 
by measuring the EW of the redshifted and blueshifted 
portions of the line separately. The wavelength regions 
adopted to sample the red and blueshifted portions of 
the line profile are listed in Table 21 and shown in Fig- 
ure [9[ where we present close-ups of the time evolution 
of each star's He I A10830 and Pa/3 line profiles. We 
selected these wavelength regions so as to best capture 
the change in morphology of the He I line on each side 
of the rest velocity while avoiding nearby contaminating 
features. 

3.3. SED analysis 

To quantify the relative contributions that the stel- 
lar photosphere and circumstellar disk make to the [3.6] 
light curve, we have constructed and interpreted broad- 
band SEDs for GY 292, 308, and 314 by combining our 
Spitzer /IRAC data with optical and NIR photometry re- 
ported in the literature. R & I band photometry is taken 
from Gordon fc Strom (1990, unp ublished; but reported 
in Table 2 of lWilking et al.ll2005f ) JHK S photometry i s 
taken from the 2MASS database (jSkrutskie et al.ll2006|) , 
and median IRAC magnitudes are adopted from the light 
curves reported here. 

The resulting SEDs for GY 292, 308, and 314 are shown 
in Figure El where we compare each YSO SED to a 0.8- 
5.0 /im diskless spectral standard of the ap propriate spec- 
tral t ype from the IRTF spectral library (Cushing et al. 
I2005D . Each standard spectrum was artificially reddened 
to simulate the extinction we mea sure toward each Y SO 
(see %T2~T1 and Table [6]) using the iFitzpatrickl (|1999t ) ex- 
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TABLE 6 
Mass Accretion Rates 



Object 


SpT 


(K) 


A H 
(mag) 


L(Br 7 ) L(Pa/3) 
(log L/Lq) 


Lacc(Br 7 ) Lace (Pa/3) 

(log L/L ) 


Afacc(Br 7 ) Mace (Pa/3) 

(IogM/M /y) 


GY 292 
GY 308 
GY 314 


K5 
MO 
M2 


4350 
3850 
3514 


2.21 
1.86 
1.12 


-4.11 -3.53 
-4.30 -3.74 
-4.89 


-0.80 -0.80 
-0.97 -1.08 
-1.50 


-8.47 -8.47 
-8.55 -8.67 
-9.19 



Computed from spectral type using the scale derived by LR99 

GY292 
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3 0.10 

E 

m 0.05 



0.00 
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?~ 

CO 
Ph 



Ratio = 0.29 
= 0.10 
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Fig. 5. — Br 7 /Pa/3 (top panels) and Pa 7 /Pa/3 (bottom panels) line ratios for GY 292 and GY 308, the two YSOs with reliably 
detected Pa/3 and Pa 7 emission. The solid line in each panel has slope equal to the derived line ratio, which we compute as 
the weighted average ratio of the da ta; la intervals a re plotted as dotted lines. Our Bry/Pa/3 line ratios are consistent with 
Br 7 /Pa/3~ 0.25, the value derived bv lBarv et all (2008) for a large sample of YSOs. This line ratio can be reproduced by models 
assuming Case B recombination theory for a wide range of temperatures and densities. There is no compelling evidence for 
variability in the either line ratio. 



tinction law with an assumed Ry = 3.1 (though the re- 
sult is largely insensitive to the adopted Ry value; see 
further discussion later in this section). We then com- 
puted integrated fluxes at JHK S and IRAC [3.6] and 
[4.5] by convolving the spectra with the appropriate fil- 
ter response curves. The standard spectrum was normal- 
ized so as to match the YSO flux at 1.24 /xm, where the 
contribution due to veiling is low. 

Assuming the spectral standards serve as reasonable 
proxies for the photospheric components of our YSO 
fluxes, the emission arising from the circumstellar disk 
at 3.6 /zm is then F3.6 — i^phot.s.ej where F3.6 is the ob- 
ject's median [3.6] flux and -F p hot,3.6 is the integrated 
flux of the reddened spectral standard in the IRAC [3.6] 
band. This modeling predicts that the disk contributes 
76% (GY 292), 59% (GY 308) and 65% (GY 314) of 
each YSO's total 3.6 /xm flux. The circumstellar disk 



is even more dominant at IRAC [4.5], where a similar 
calculation implies disk fractions of 91% (GY 292) and 
85% (GY 314; note that we did not obtain [4.5] data 
for GY 308). We also estimated the expected amount 
of veiling (r x = F disktX / F pho t,\) in the H and K s bands 
by comparing the filter-integrated fluxes of the reddened 
spectral standards to the YSO 2MASS fluxes, deriving 
{r H ,r K } veilings of {0.4,1.3}, {0.2,0.8}, and {0.4,0.7} 
for GY 292, GY 308, and GY 314, respectively These 
values should be considered lower limits, as they assume 
negligible veiling at J band (discussed further later in 
this section), and they are consistent (as lower limits) 
with the more robust results we derived through direct 
measurements of absorption lines in ^3.2.21 

These SED-based veiling/disk emission estimates are 
subject to three dominant sources of uncertainty: the 
adopted extinction law, potential presence of J band veil- 
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Fig. 6.— Optical/IR SEDs for GY 292, 308 and 314 (di- 
amonds), including RIJHK S from the literature and our 
Spitzer IRAC photometry. Each SED is compared to the 0.8- 
5.0 /im spectrum of a diskless fie ld star of the sam e spectral 
type, artifically reddened with the lFitzp atrick 1999 extinction 
law (with Rv ~ 3.1) to match the target star's spectroscop- 
ically inferred Ah, and normalized to the match the target 
star's J band flux. GY 292, 308, and 314 all possess near- and 
mid-infrared excesses relative to their reddened spectral stan- 
dards, implying that about 76%, 59%, and 64% of their [3.6] 
emission, respectively, arises from their circumstellar disks. 

ing, and long-term photometric variability. However, as 
we demonstrate below, these effects do not appear large 
enough to alter our basic conclusion that, for each each 
YSO, the mid-IR flux is dominated by disk emission: 

• Adopted extinction law: Our SED analysis depends 
centrally on the artificially reddened spectral stan- 
dard that serves as a proxy for flux expected from 
the YSO's stellar photosphere. By directly affect- 
ing the shape of that assumed photospheric flux, 
errors in either the magnitude or wavelength de- 
pendence of the assumed extinction will translate 
directly into errors in our inferred veilings or disk 
emission fractions. Several previous studies, how- 
ever, have established that the NIR extinction law 
is highly insensitive to dust pro perties, or equiva- 
lently to the choic e of Ry (e.g.. iFitzpatrickl 1999; 
llndebetouw et alJl2005HChapman et al.ll2009l) . We 
confirmed this result by re-calculating disk fraction 
and veiling estimates after reddening the standard 
with an Ry = 5.5 (i.e., grayer) extinction law, and 
the resulting [3.6] disk fractions decreased by only 
1-2%, while veiling estimates decreased by no more 
than Ar\ = 0.1. 

• Potential J band veiling: Our SED analysis im- 
plicitly assumes that each YSOs J band flux is 
entirely photospheric, i.e., rj = 0. The estimate 
we infer for the fraction of each YSO's flux that is 
non-photospheric in nature is therefore quite con- 
servative: if a YSO does indeed have a nonzero 
J band veilin g (as most do: m any CTTSs possess 
Tj ^0.1-0.2; ICieza et al.l 120051 ). the disk excesses 
inferred from our SED based analysis will be un- 
derestimated by a factor proportional to the un- 



acknowledged J band excess. Specifically, to cor- 
rect our estimates for a YSO with a typical J band 
veiling of 0.2 (corresponding to a flux correction 
of about 15-20%), we would need to increase tr 
by 0.2-0.3, r K by 0.3-0.5, and the [3.6] disk frac- 
tion estimate by 4-8%. Thus, the results we derive 
above should be considered robust lower limits on 
disk fraction and veiling, and the mid-IR emission 
from the circumstellar disk is likely to be even more 
dominant at 3.6 /im than we predict in our sample 
of YSOs. 

• Long term photometric variability: The accuracy of 
the broadband SEDs we have constructed is lim- 
ited by the amplitude of each YSO's long-term 
photometric variability, due to the approximately 
ten-year timespan between the 2MASS JHK S and 
our Spitzer IRAC measurements. Such long-term 
variations are indeed observed, thoug h not ubiq- 
uitous, in young stars; for example, iParks et al.l 
(|in prep. I ) find that 31% of the variable stars in 
their study of p Oph exhibit modest to significant 
K s band variations over timescales of a few months 
to 2.5 years, although the long-term variability am- 
plitude is limited to < 0-3 mag in the majority of 
cases. Furthermore, IScholzl ([20121 ) find that there 
are very few YSOs (2-3% of cluster members in sev- 
eral nearby star-forming regions including p Oph) 
that are highly variable (> 0.5 mag) over years- long 
timescales. To test the potential effect of long-term 
variability in our sample, we compared the 2MASS 
magnitudes of our objects (listed in Tabled]), which 
were measured between 1997-2001 , to the pho- 
tomet ry obtained in 1993-1994 by iBarsonv et al.l 
(|1997f ). and find changes in all thr ee NIR bands 
of only about 0.1 mag on average. IScholzl ((20]1) 
predict that the absolute error in a particular pho- 
tometric measurement of a YSO due to long-term 
variability is 5-20%, with some dependence on age 
and the fraction of emission produced by the disk. 
As the discussion of potential uncertainties due to 
J band veiling demonstrates, errors on the 15-20% 
level limit the precision of the exact disk fraction 
values derived above, but the qualitative result that 
the majority of [3.6] emission originates in the cir- 
cumstellar disk is robust against potential long- 
term photometric variations. 

4. DISCUSSION 

We now seek to use the time-averaged and time- 
resolved spectral diagnostics measured in the previous 
section (e.g, H & K s band veiling, Br7, Pa/3, and 
He I A10830 EWs, and spectrally integrated J — K s col- 
ors) to inform our understanding of GY 292, 308, and 
314's mid-IR light curves and evolutionary states. After 
a brief discussion of the implications of the time-averaged 
properties of our YSOs, we test for signatures of the phys- 
ical mechanisms driving the moderate mid-IR variability 
in our sample by computing the linear Pearson's correla- 
tion coefficient between their IRAC light curves and the 
relevant spectral diagnostics. The time dependence of 
our full suite of spectral and photometric measurements 
is summarized for these three YSOs in Figure [7J 
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Fig. 7. — Evolution of all time-dependent quantities mea- 
sured from the light curves and spectra of GY 292, 308, and 
314: IRAC 1 (i.e. [3.6]) photometry (black), IRAC (1 - 2) 
color (red), K band veiling (blue), Br7 EW (accretion rate 
diagnostic, orange), and spectrally integrated J — K„ color 
(extinction diagnostic, purple). The average value is listed to 
the right of each time series, and indicated on the plot by a 
dashed line; the plot displays fractional deviations from this 
average. The magnitude of variations in tk and Br7 EW have 
been suppressed by a factor of two for clarity. GY 308 and 
314 exhibit significant mid-IR variability, and GY 292 and 
308 exhibit variability in the accretion-sensitive Br 7 emis- 
sion lines; we identify no significant correlations between the 
time dependence of these physical diagnostics, however. 



4.1. Evolutionary State 

Theory and observations suggest that YSOs begin their 
lives as strong accretors with massive disks, with the ac- 
cretion rate and disk mass declining over time as cir- 
cumstellar material is depleted by accretion, dispersed 
by stellar outflows, or collected into dense planetesi- 
mals. All three of GY 292, GY 308, and GY 314 
we re classified as Class I I sources (i.e. T Tauri stars) 
bv iBontemps et all (|2001l k although we lack the longer 
wavelength mid-IR coverage to reliably compute a, the 
spectral index commonly used to classify YSOs, the 
SEDs constructed and disk fractions computed in £13.31 
are qualitatively consistent with this classification (see 
also Figure |6|). Star formation is a dynamic process, with 
notable examples of young stars whose M acc and veiling 
have undergone significant chan ges (e.g. iGreene &; Ladal 
H99tl IConnellev fc Greenel l2lH0l) . We therefore define a 
young star's "evolutionary state" as a measure of the 



source's characteristic levels of veiling and mass accre- 
tion, such that heavy accretors/ veiled young stars are 
less "evolved" than their more quiescent kin. This evolu- 
tionary sequence will not, in general, map well to a star's 
chronological age, but instead to its present day levels of 
disk accretion and emission. The time-averaged veiling 
and mass accretion rates of our three YSOs differ mea- 
surably, such that our data suggest a clear "evolution- 
ary sequence" that coincidentally follows their numerical 
designation: GY 292 has the highest accretion rate and 
veiling, GY 314 has the lowest accretion rate and veiling, 
and GY 308 is intermediate between the two. 

Veiling has been shown to correlate strongly with the 
evolutionary state of a forming star: the most embedded 
objects are shrouded in so much circumstellar material 
that photospheric features can disappear entirely, while 
disk emission decreases gra dually as matter is ac creted 
onto the evolving star (e.g.. IGreene fc Ladal ll996). Our 
veiling measurements displayed graphically in Figure HI 
agree well in both scale and spectral shape with the disk 
excesses implied by the SED analysis in ^3.31 and identify 
GY 292 as the most heavily veiled (and thus the YSO 
having the most dominant circumstellar disk) of these 
three; GY 314 is the least veiled, but only modestly less 
so than GY 308. 

The evolutionary sequence suggested by the level of 
veiling in our YSO targets is also consistent with the 
mass accretion rates measured in Section r3.2.3l and shown 
in Table El GY 292 has the highest time-averaged mass 
accretion rate (10 -85 M yr _1 ), and GY 314 has the 
lowest (10 -9 2 M Q yr _1 ). Additionally, the amplitude of 
variability in M acc (as inferred from the standard devi- 
ation of the EW(Br7) time series, <7Br m Figure [7]) is 
largest in GY 292 and smallest in GY 314, hinting that 
higher characteristic accretion rates may also be accom- 
panied by more pronounced short-timescale M acc varia- 
tions. We do note, however, that while GY 292's near- 
infrared spectrum exhibits notable variability, particu- 
larly in the H I and He I accretion diagnostics, it is rela- 
tively quiescent in both IRAC bands over the Spitzer ob- 
servation window. This suggests that, despite the broad 
evolut ionary trends identified by Moralcs-Caldcro n et al.1 
(| 2 1 It) . mid-IR variability is not a uniform, monotonic 
function of accretion or disk activity, at least on week- 
to month- timescales. 

The mass accretion rate (as inferred from BT7 EW) of 
each YSO in our sample varies by a factor of 3-5 over our 
observation window, while time-averaged accretion rates 
differ by about a factor of five across the sample. Since 
the latter are averaged quantities, this suggests a real 
difference in evolutionary state of the three YSOs that 
is more significant than the short-term, stochastic varia- 
tions observed in any one YSO. While this is not strict 
proof of evolutionary ordering - indeed, Class II YSOs 
have been seen to exhibit differences in EW(Br7) and 
EW(Pa/ 3) of up to a factor o f ~10 across a flux-limited 
sample ([Greene fc Ladal[l996l ) - the fact that our objects 
seem to be measurably more disparate in their character- 
istic accretion activity than any single YSO is variable, 
combined with a clear difference in the level of veiling 
across the sample, leads us to propose the evolutionary 
sequence above as the simplest (but not only) explana- 
tion that is consistent with our measurements. 
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4.2. Possible explanations for mid-IR variability 

In the following sections, we explore if each YSO's ob- 
served IRAC [3.6] variability can be attributed to phys- 
ical processes that correlate with the its veiling, extinc- 
tion, and/or accretion rates. 

4.2.1. Veiling correlation 

Detailed radiative transfer models indicate that a 
YSO's NIR and mid-IR excess emission both likely arise 
from the inner most regions of its circumstellar disk. 
lEspaillatl (|2009[ ) showed that more than 90% of the 7 /mi 
flux produced by a typical circumstellar disk originates 
inward of ~0.25 AU; excess emission at shorter wave- 
lengths will be even more centrally condensed, arising 
from the warmest material at the smallest radii. As disks 
within systems with M acc of 10 -9 to lO~ 8 M0/yr are typ- 
ically i Jxuncatec^ 0.15 AU 
(e.g.. lEspaillat et al.ll2007HLe Blanc et al.ll2011l) . this im- 
plies that the vast majority of both the K band and 
3.6 /Ltm flux is emitted within the inner ~0. 05-0. 20 AU 
of the disk, likely from the inner wall at the disk's trunca- 
tion radius. Thus one might naturally expect [3.6] mag- 
nitude and NIR veiling to be correlated. 

To test this possibility in our sample, we explicitly look 
for a correlation in the time evolution of each YSO's [3.6] 
flux and NIR veiling rjj and rjc- Interpolating the IRAC 
[3.6] light curve onto the spectroscopic epochs, we cal- 
culate the linear Pearson's correlation coefficient, r, be- 
tween the two time series; we find r = 0.59 and 0.38 (GY 
292), 0.30 and 0.21 (GY 308), and 0.31 and 0.23 (GY 
314), where the two rvalues listed for each object are for 
ru vs. [3.6] and tk vs. [3.6], respectively. These values 
are well below the 95% confidence value of |r| = 0.811, 
and thus we are unable to confidently link the observed 
mid-IR variability with changes in H and K band veiling. 
However, veiling variability is only present at a modest 
level in all three objects (in contrast with accretion vari- 
ability, which is clearly evident in at least two YSOs; see 
g4.2.3|) . Thus modest changes in disk emission could still 
potentially explain the IRAC light curves of marginally 
variable objects such as GY 292. 

4.2.2. Extinction model 

We next investigate whether changes in each YSO's 
line-of-sight extinction could explain the observed mid- 
IR variability. To explore this possibility, we construct a 
simple extinction model with three assumptions: (1) the 
median IRAC [3.6] magnitude represents the source's me- 
dian state; i.e., the spectral parameters inferred earlier 
and published 2MASS magnitudes are a good represen- 
tation of that epoch; (2) the variability in the measured 
IRAC [3.6] magnitude is due solely to changes in extinc- 
tion; and (3) the inters tellar extinction law derived by 
llndebetouw et al.l (|2005l ) is an appropriate proxy for disk 
extinction^- We parameterize extinction as J — K s color, 
and compare the model predictions to the spectrally in- 
tegrated J — K s colors computed in £13.2.11 

In particular, we used our [3.6] measurements with 
the above assumptions to calculate the extinction-based 

5 The Indcbctouw ct al. (2005) and Fitzpatrick (1999) extinction 
laws agree to within a few percent in the infrared, thus we adopt 
the former for this model since it is native to the IRAC filter set. 



model Aj and Ak s at each IRAC epoch, then com- 
puted the J — K s color predicted by the model as 
(J - K s ) modcl = (J - IQ2MASS + {Aj - A Ks ). We 
compare these model colors to the spectrally integrated 
J — K s colors computed in i j3.2.1l in Figure |SJ As the fig- 
ure demonstrates, the spectrally integrated J — K s colors 
show no significant variability across the spectroscopic 
epochs, limiting extinction-driven changes in J—K s color 
to less than ~0.2 mag for all three YSOs. Modeling the 
[3.6] light curves as fully due to extinction variations, by 
contrast, predicts J — K s variations at the 0.2-0.8 mag 
level over the full duration of the photometric monitor- 
ing. The J — K s variations predicted by the extinction 
model are definitively inconsistent with GY 314's spec- 
trally integrated J — K s colors, ruling out extinction as a 
potential mechanism to explain the structure observed 
in GY 314's [3.6] light curve. GY 292 is not signif- 
icantly variable at [3.6], by contrast, constraining any 
potential extinction variations to < 0.2 mag in J — K s 
for this YSO. GY 308 is an intermediate case; while the 
model and spectrally integrated J — K s are largely con- 
sistent at most of the spectral epochs, there are spectro- 
scopic datapoints that disagree with the predictions of 
the extinction-based model at the 1.5 a level, suggesting 
that mechanisms other than extinction must be at least 
partially responsible for the observed mid-IR variability. 

The weakest link in the above calculation is the as- 
sumption that an interstellar extinction law accurately 
replicates the extinction due to dust in a circumstel- 
lar disk. Extinction curves for circumstellar dust are 
often thought to be gray er than for interstellar dust 
(e.g., iFlahertv et all [20071) . so a model based on an in- 
terstellar extinction law could potentially overestimate 
the scale of the expected NIR extinction variations re- 
quired to reproduce a given [3.6] light curve due to ex- 
tinction from intermediate disk structures. Nevertheless, 
the exact wavelength dependence of extinction in disks 
is not well-constrained and dust properties could vary 
significantly with e volutionary age, s tellar properties, or 
viewing angle (e.g.. iSung et al.l l2009). so the above cal- 
culation serves as a well-defined test of an extreme case 
of disk extinction as the source of mid-IR variability. 

4.2.3. Accretion rate correlation 

Enhanced accretion activity could potentially heat the 
inner disk, both through viscous dissipation and via en- 
ergetic radiation emitted by the accretion shock, and 
thereby produce an increase in the mid-IR emission from 
the disk's inner edge. Physical changes in the disk struc- 
ture related to the accretion process could also alter the 
disk's emission properties. If active, these mechanisms 
should produce a correlation between the YSO's [3.6] 
magnitude and the strength of accretion-sensitive lines, 
such as Br7. 

As discussed in i)4.11 the mass accretion rate, as traced 
by Bi'7 EW, appears to vary significantly with time in our 
sample, especially in GY 292 and GY 308. To formally 
quantify that EW(Br7) is variable in our YSOs, we adopt 
the null hypothesis that Br7 EW is constant in time, 
and perform a % 2 minimization to fit our measured B1-7 
EWs with a constant accretion rate model. This model 
produces reduced x 2 values of 29, 28, and 3.7 for GY 292, 
308, and 314, respectively. The significant deviation from 
the time-independent accretion model is strong evidence 
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Fig. 8. — Evolution of J — K a color (derived by integrating 
the SpeX spectra over the 2MASS filter bands; large red cir- 
cles) superposed on the variations predicted for the J — K s 
color from a simple model in which variations in each YSO's 
[3.6] light curve are solely attributed to changes in the YSO's 
extinction, Ay 3 ^ (small black circles). Extinction variations 
could induce GY 292's low level of [3.6] variability, but at- 
tributing GY 314's observed [3.6] variations to variable ex- 
tinction predicts variations in J — K s colors that are incon- 
sistent with the observed SpeX data. Modeling GY 308's 
[3.6] varations as extinction events predicts J — K a variations 
that are only marginally consistent with the observed SpeX 
spectral data. 

that these YSOs undergo measurable changes in M acc on 
day- and week- timescales, particularly the presumably 
less evolved YSOs GY 292 and 308. 

To test if the variability in the accretion rate is re- 
lated to the structure in the [3.6] light curve, we calcu- 
late the level of correlation between each YSO's observed 
Br7 EW and [3.6] flux, following the procedure of §4.2.11 
above. We calculate rvalues of r = -0.31 (GY 292), 
-0.59 (GY 308), and -0.41 (GY 314), each well below 
the threshold values |r| = 0.811 (6 data points) or 0.878 
(5 data points) that indicate 95% confidence in a correla- 
tion. We are thus unable to establish that the observed 
changes in mass accretion are primary contributors to 
each YSO's mid-IR variability. 

4.3. #e/A10830 

We conclude this study by examining the temporal cor- 
relation between the Pa/3 emission line and the structure 
of the He I A10830 line. Previous observational and the- 
oretical work has shown that YSO He I and H I line 
profiles can both be shaped by accreting and outflow- 



ing material (e.g.. lAlencar et ai]|2005l iKwan et al. l l2007l; 
Kurosawa et alj|201lD : nonetheless, the typical emergent 
profiles seen in these lines from YSOs differ significantly, 
as do their assumed sensitivity to accretion and outflow. 

Pa/3 is seen in emission in most Class II YSOs, with a 
minority exhibiting a mod est amount of sub-conti nuum 
redshifted absorption (e.g., Folha & Emerson 2001 found 
that 42/50 CTTS exhibited Pa/3 in emission, of which 
13 also exhibited redshifted absorption). Pa/3 line 
strengths also correlate quite well with mass accretion 
rates inferred from blu e continuum excesses (e.g., Fig. 3, 
iMuzerolle et al.lll998| ). such that it is often adopted as 
an accretion indicator, despite the fact that much of the 
observed emission could plausibly arise from outflowin g 
material (see, e.g., Fig. 9. lKurosawa k, Romanoval f2012). 

He I A10830, in contrast to Pa/3, exhibits absorp- 
tion components just as frequently as emission: in their 
analysis of high resolu t ion H e I A10830 profiles of 39 
CTTSs, lEdwards et ail ((2006h found 35 and 33, respec- 
tively, that exhibited absorption and emission compo- 
nents. Studies of He I A10830's excitation conditions 
indicate that outflows are likely responsible for most of 
the emission, with outflows and accretion responsible 
for blueshifted and redshifted absor ption , respectively 
(|Fischer et al.ll2008t IKwan fc Fischerll201lD . It has been 
suggested, however, that the outflows t raced by He I 
may themselves be accretio n-driven (e.g., E dwards et al.l 
120061 : IMatt fc P udritz 2008), suggesting that correlations 
may exist between the (presumably) accretion-sensitive 
H I lines, such as Pa/3, and the predominantly outflow- 
sensitive lines, such as He I A10830. 

To test for potential correlations between the kincmat- 
ically simple Pa/3 emission line, and the kinematically 
richer He I A10830 line, we calculated linear Pearson's 
correlation coefficients (r) between the EW of the full 
Pa/3 line and the EWs of the blueshifted and redshifted 
components of He I (computed in the regions bounded 
by the black dotted line, which is at the systemic rest ve- 
locity, and the blue and red dotted lines, respectively, in 
Figure [9]) . These constant- width windows were selected 
such that the entire P-Cygni profile, where present, was 
captured at all epochs. Pa/3 and He I EWs were corrected 
for the effects of veiling using the procedure of ^3.2.2l and 
measurements of the Al I AA13123, 13151 doublet in both 
our target and standard star spectra. The results of these 
correlation calculations are shown visually in Figure 1101 

We find a clear correlation between the Pa/3 EW and 
the EW of the red side of He I A10830 for all three YSOs 
(r = 0.79, 0.79, and 0.74 for GY 292, 308, and 314, 
respectively), but only a modest (and not statistically 
significant) correlation between the blueshifted side of 
this He I line and Pa/3 emission (r = —0.59, 0.43, and 
0.28; |r| > 0.6 indicates statistical significance for this 
sample size). For all three YSOs, Pa/3 emission increases 
more or less monotonically with the strength of the red- 
shifted component of the He I A10830 line, suggesting 
that both of these line components are indeed sensitive 
to accretion activity, consistent with the interpretation 
advanced bv lFischer et al. (2008f). This could reflect ei- 
ther a) the presence of increased He I A10830 emission 
from winds at epochs with stronger Pa/3 emission, which 
itself could be due to a combination of emission from ac- 
cretion or outflows, or b) a geometrical effect whereby 
emission from Pa/3 is enhanced during epochs when the 
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Fig. 9. — He I A10830 and Pa/3 profiles for each of our targets. The blueward (redward) He I EW was computed over the 
spectral region between the dotted blue (red) line and the dotted black line at the systemic rest velocity. All spectra were shifted 
so as to center the Pa/3 line center at its rest wavelength. Redshifted He I emission appears to be correlated with Pa/3, while 
blueshifted absorption is less so (see also Figure [10)) . The epoch is MJD— 55000. 
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Fig. 10. — Linear fits to the relations between Pa/3 equivalent width (an accretion diagnostic) and each of He I A10830 
blueshifted absorption (left panels) and He I A10830 redshifted emission (right panels) for the YSOs GY 292, GY 308, and GY 
314. Correlation coefficients r with magnitude greater than about 0.6 indicate statistically significant correlations, suggesting 
that the redshifted portion of He I may be more closely correlated with accretion activity in these YSOs than the blueshifted 
portion. Note that negative EWs indicate emission, while positive EWs measure absorption. 



magnetospheric accretion column is perpendicular to our 
line of sight, and therefore is also unable to produce sig- 
nificant redshifted absorption in He I A10830. 

We note, however, that as each of our targets is a "typi- 
cal" Class II YSO, the correlation we detect between Pa/3 
emission and redshifted He I A10830 may only hold for 
sources over a limited range of evolutionary states, or 
with a particular source geometry. Indeed, FU Ori stars 
demonstrate a clear case where this correlation cannot 
hold, for example, as these sources lack the H I emission 
that serves as a proxy for accretion rate for the YSOs 
in our sample (|Connellev & Greene! [20Toh . The FU Ori 
phenomenon is sufficiently distinct from the "standard" 
magnetospheric accretion mechanism thought to operate 
in our Class II sources, however, that the absence of this 
correlation in the FU Ori phase is not dispositive of the 
implied connection between the Pa/3 emission, redshifted 
He I A10830 absorption, and the underlying accretion 
mechanism. 

The fact that the blueshifted component of He I is not 
significantly correlated with accretion activity over time 
seems, at first blush, to argue against a model wherein 
YSO winds are instantaneously accretion-powered. How- 



ever, despite the lack of a formal statistical correlation 
between the Pa/3 and blueshifted He I EWs, visual in- 
spection of Figure [TU] does reveal a qualitative (although 
modest) connection between these quantities. Addition- 
ally, when considering the three YSOs as a sample, the 
time-averaged strength of the He I A10830 blueshifted 
absorption feature is correlated with the time-averaged 
accretion rate as inferred from Pa/3 and Br7 emission, in 
the sense that GY 292 has both the highest M acc and 
the largest time-averaged blueshifted He I EW, while 
GY 314 has the lowest in both. Furthermore, we do find 
a statistically significant correlation between the total 
He I A10830 line strength, defined as |EW rcd | + |EW b i ue |, 
and the Pa/3 EW in these YSOs, deriving r = —0.69 
and r = -0.65 for GY 292 and GY 308, respectively (we 
again excluded GY 314 because it lacks Pa/3 emission at 
se veral epochs). This re sult agrees with the trend found 
by Edwards et all (|2006| ). who noted an increase in both 
total He I A 10830 line strength and Pa7 emission with 
increasing Y band veiling (and thus a mutually increas- 
ing relation between He I and Pa^0). Interestingly, the 

6 We measure Pa/3 instead of Pa7 due to our greatly increased 
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slope of the relation between blueshifted He I A10830 
absorption and Pa/3 emission reverses sign (as indicated 
by the sign of r) for GY 292 as compared to GY 308 
and GY 314, which suggests further complications in the 
dynamics of the accretion/outflow relationship. This dif- 
ference in sign may represent further differences in accre- 
tion or outflow geometry between these YSOs, or since 
the correlations are not statistically significant, the re- 
versal in sign may simply be due to statistical noise. 

Detailed time-series studies of more objects, and at 
higher spectral and temporal resolution, are needed to 
more fully explore the apparently complicated relation- 
ship between stellar outflows and accretion activity. Di- 
rect, simultaneous comparisons of the He I A10830 pro- 
file with other accretion and outflow indicators, such as 
UV-continuum excess, could also help clarify the relative 
contributions of accretion and outflow to this complex 
line profile. 

5. SUMMARY 

1. We derived spectral types, extinctions, veilings, 
and mass accretion rates for a small sample of 
YSOs in p Oph. All three are Class II YSOs 
(CTTSs) and active accretors with time-averaged 
mass accretion rates of ~ 10 -9 2 — 10 -85 M Q /yr. 
Our results suggest an evolutionary sequence 
(which may not directly map to an age sequence) 
from least to most evolved of GY 292, GY 308, 
then GY 314. 

2. IRAC mid-IR light curves were also obtained for 
these objects as well as the Class I YSOs YLW 15a 
and YLW 16b. Three of five of these YSOs ex- 
hibit statistically significant variability with [3.6] 
amplitudes of > 0.2 mag. Interestingly, GY 292, 
the Class II YSO that is the least evolved (i.e. the 
strongest accretor and most veiled) is the least vari- 
able in the mid-IR over day- to week- timescales, 
suggesting that mid-IR variability amplitude is not 
monotonic as a function of accretion activity or 
dominance of the circumstellar disk. 

3. To investigate possible mechanisms connected with 
mid-IR variability, we examined the extinction, 
veiling, and mass accretion rate of each YSO as 
a function of time. Accretion rates, as traced by 
H I emission, are seen to vary by as much as 
a factor of several over our two month baseline, 
while variations in the emission by the inner disk, 
as traced by NIR veiling, are much more mod- 
est. However, we cannot directly link changes in 
either accretion or disk emission to variations in 
the IRAC light curves. Extinction is relatively con- 
stant over our observational baseline, and we can 
rule it out as the sole cause of mid-IR variability 
in at least one object, GY 314. Self-obscuration by 
the disk remains a possible scenario in the other 
two YSOs; this possibility could be further investi- 
gated through higher precision mid-IR photometry 
and/or polarimetry to detect scattered light 

S/N at the wavelength of the former as compared to the latter; the 
relative strengths of these t wo rai o wcr = 3 transitions appears to be 
consistent in our data (see 813.2.3)1 . 



4. We find that redshifted He I A10830 emission, 
where present in our spectra, shows both quantita- 
tive and qualitative similarities to the accretion- 
sensitive H I emission lines, suggesting that the 
winds presumably traced by the He I line may be 
connected with accretion activity. Blueshifted He I 
absorption, while not formally correlated with Pa/3, 
is strongest (in a time- averaged sense) in GY 292, 
the least evolved source, and practically absent in 
GY 314, the most evolved source. Additionally, the 
total He I A10830 line strength is correlated with 
the Pa/3 EW, further supporting a connection (al- 
beit a complicated one) between accretion, veiling, 
and winds in YSOs. 
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APPENDIX 
ADDITIONAL SPECTRA 

We present the SpeX K band spectra of our secondary targets, YLW 15a and 16b, in Figure [TT] 
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FlG. 11. — K band spectra for the two secondary targets in our sample. Both objects have rising K band slopes, and so the continuum 
was fit with a second-order polynomial. 
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